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Abstract 



We study conformal field theory correlation functions relevant for string dia- 
grams with open strings that stretch between several parallel branes of different 
dimensions. In the framework of conformal field theory, they involve boundary 
condition changing twist fields which intertwine between Neumann and Dirichlet 
conditions. A Knizhnik-Zamolodchikov-like differential equation for correlators 
of such boundary twist fields and ordinary string vertex operators is derived, and 
explicit integral formulas for its solutions are provided. 
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1. Introduction 



D-branes [30] have become the most important ingredient of the new picture of string 
theory that has emerged in recent years. They have shaped a new understanding of non- 
perturbative effects in string theory and of low-energy effective theories associated with 
string theories. In the latter context, systems of many branes are of particular importance, 
since they provide a natural way to include non-abelian gauge theories into string theory 
[40]. Systems of several branes of different dimensions, most notably of Dl- and D5-branes, 
play a major role in proposals of how to derive the Bekenstein-Hawking entropy of black 
holes from string theory [39,27]. More recently, stacks of branes and anti-branes have been 
reconsidered in connection with a i^-theoretic classification of branes [41] , based on results 
concerning tachyon condensation in [37]. 

Some qualitative features of such systems can be uncovered within a target space approach. 
But e.g. the process of brane-antibrane annihilation in the last-mentioned application, or 
the properties of near-extremal black holes involve the analysis of non-BPS states for 
which the world-sheet approach is better suited, since it does not critically depend on 
supersymmetry. Computation of CFT correlation functions is indispensable if one wants 
to deal with problems like Hawking radiation off Dl-D5-systems, or for a clean discussion 
of bound state formation [23] . 

The world-sheet formulation of string sectors with branes is well known in string theory, 
although mainly in connection with flat targets. The general setup involves boundary con- 
formal field theory as introduced and developed by Cardy [8,9,10] and first exploited for 
string theory by Sagnotti [34] . CFT on surfaces with boundaries exhibits a very interesting 
internal structure and finds interesting applications beyond string theory: For many s-wave 
dominated scattering processes, the universal behaviour is described by a boundary CFT 
in two dimensions, irrespective of the dimensionality of the original system. The most 
famous problem that could be tackled with boundary CFT methods is the Kondo effect in 
condensed matter physics [1]. 

In string theory, methods of boundary CFT are not only valuable in the study of situations 
without the BPS-property, but also to uncover non-classical features like unexpected mod- 
uli [37,33] and non-commutative geometry; see e.g. [15,35,36,3,16] and references therein. 
Moreover, they allow one to analyze D-branes in non-geometric string compactifications 
such as Gepner models [29,38,32,24,6,13]. 

In this paper, we ask how to compute CFT correlators describing string amplitudes of 
arbitrary closed and open string vertex operators in the presence of multiple flat branes 
in JBL D . The open strings involved stretch between two or more different branes, which 
may have different dimensions. It appears that no systematic method for the computation 
of those string diagrams, which contribute to scattering processes in higher orders of the 
string coupling constant, is available in the literature. See, however, [25] for some sample 
computations of scattering amplitudes in the presence of a pair of branes. 

The world-sheet description requires surfaces with several boundary components. We 
restrict our attention to diagrams without internal closed string loops so that we can 
map the world-sheet to the disk or to the upper half-plane, but with different boundary 
conditions assigned to consecutive intervals on the boundary; see Figure 1. We focus on 
parallel branes here; thus we can reduce our analysis to a one-dimensional target. Results 
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for R D with D > 1 follow by taking tensor products. The boundary state for a p-brane 
involves p + 1 Neumann and D — p — 1 Dirichlet boundary states of single free bosons. 
The interesting transitions between boundary conditions in a one-dimensional target are 
those from Neumann to Dirichlet or vice versa. They are mediated by boundary fields of 
a special type, namely boundary condition changing twist Gelds. 




Figure 1: The upper half-plane with a sequence of Neumann (solid intervals) 
and Dirichlet (dashed intervals) boundary conditions along the real line. The dots 
between Neumann and Dirichlet intervals mark insertions of boundary condition 
changing twist fields, while crosses on the boundary or in the interior refer to 
insertions of ordinary open resp. closed string vertex operators. Such a world- 
sheet diagram can be understood as Hawking radiation (closed string states) from 
a system of branes (multiple changes of boundary conditions) with simultaneous 
inner excitations (open string states). 

Conformal boundary conditions which preserve the chiral algebra W of the theory are 
parametrized by certain automorphisms O of W, together with the amplitudes of one- 
point functions [32]. If W is the U(l) current algebra, O can act as ±id on the currents, 
and the 1-point functions determine the location of a brane. If the boundary condition is 
constant along the boundary, arbitrary n-point functions can be expressed in terms of the 
usual conformal blocks of W; see e.g. [8,32,19]. 

If the gluing conditions described by an automorphism Q of W remain constant along the 
boundary the computation of correlation functions is not, in principle, a difficult problem. 
Otherwise, the simple Ward identities for the symmetry algebra W are broken, and one 
has to find new methods to construct the "twisted chiral blocks" involving a new type of 
boundary condition changing operators which correspond to twisted rather than ordinary 
representations of W. It is the aim of this article to develop a convenient formalism for 
computing such correlation functions in the case of a flat target space. 

The plan of this paper is as follows: In the next section, we look at correlation functions 
which contain just one insertion of a boundary twist field. We shall provide a complete 
operator construction of the boundary CFT, from which one can derive correlators with 
an arbitrary number of closed string vertex operators inserted in the bulk. When there are 
more than two twist fields on the boundary, such techniques are no longer available. Our 
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strategy is then to derive Ward identities for the correlation functions. They will lead us 
to Knizhnik-Zamolodchikov-like differential equations which describe the effect of moving 
insertion points for bulk and boundary fields in terms of a flat connection. We explain this 
idea in Section 3 and exploit it in the fourth section to give explicit integral formulas for the 
correlators. While some of the technical steps in setting up the Knizhnik-Zamolodchikov 
equations are rather involved, parts of our final results can be related to electrostatics. 
Section 5 comments on possible generalizations and applications. 



2. Operator formalism for a single twist field insertion 

As a simple example, we consider open strings propagating freely in the target R, with 
Dirichlet boundary conditions imposed at one end of the string and Neumann boundary 
conditions at the other. We thus have to deal with a free bosonic field X(t, a) defined for 
space variables a G [0, n] and subject to 

d t X(t,0) = , d a X(t,ir) = Q 

for all t G R. Mapping the strip to the upper half-plane H by z = exp(t + ia), X(z,z) 
satisfies Dirichlet boundary conditions for z G R>o and Neumann boundary conditions for 

z G R <0 : 

X(z,z) = x for z = z>0 , (d-B)X(z,z) = for z = z < . (1) 

Our task is to compute correlation functions in this bosonic theory which involves two 
insertions of twist fields on the boundary (see [12] for an early treatment of that problem). 
Conformal symmetry allows us to place these boundary condition changing operators at 
x\ = and x-i = oo. 

We propose to construct operators X(z, z) satisfying (1), as well as open and closed string 
vertex operators, and then to derive differential equations on the correlation functions from 
the algebraic properties of these operators and from the symmetries of the theory. First, 
we have to determine the space our fields are to act on. 

2.1 The spectrum of boundary twist fields. The space in question is spanned 
by excited states of open strings stretching between a Neumann and a Dirichlet boundary 
condition - hence the name "boundary condition changing operators" for the boundary 
fields uniquely associated to these states. There is a relatively simple technique to de- 
termine the spectrum of boundary condition changing operators that intertwine between 
two constant conformal boundary conditions B\ = (Qi,a\) and B 2 = (^2,012) in some 
boundary CFT. The state space of this boundary theory is denoted by Hi2- Because of 
the state-field correspondence, the spectrum of boundary condition changing operators is 
described through the partition function of the boundary theory, 

Z 12 (q) = Tr Hl2 q H(H) where r = L - ± . 
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By an interchange of space and time coordinates ( "world-sheet duality" ) , the open string 
1-loop diagram underlying Z\i may be viewed as a closed string tree diagram, i.e. 

Z 12 (q) = {B^qt*™ \B 2 ) where = L< P) + L ^ - ± , 

where q = exp(— 2-Ki/r) is related to the variable q = exp(27rir) as usual. The closed 
strings propagate between the boundary states \Bi) = associated with the boundary 

conditions Bi. They allow to transfer boundary conditions from the upper half-plane 
(where the Hamiltonian is H^ H) ) into a CFT on the full plane (with Hamiltonian -ff (p) ), 
see [26,9]. 

The boundary states implementing Dirichlet and Neumann conditions along the whole 
boundary are of course well known, see e.g. [7]. Let a^ p) , a^ p) be two commuting sets of 
oscillator modes (in the plane CFT) with standard commutation relations. The ground 
states | k) of their Fock spaces are labeled by the momentum fcel. Neumann and Dirichlet 
boundary states are given by 

|iV) = i=exp{-^ia^a^}|0), 
V2 ^ n 

, ^ ! (2) 

| £>(xo) )= dk exp{ £ - a™ a™ } \k) , 
J „>i n 

where x G M, as in (1), denotes the location of the "D-brane", i.e. x |_D(x )) = x |-D(x )) 
for the center of mass coordinate x. 

If we build the boundary state \p) for a p-brane in M. D as a tensor product of \N) 
(j> + 1 times) and \D(xq)) (i = p + 2, . . . , D) from eq. (2), the partition function Z w (q) 
counts boundary fields that do not change the boundary condition. They describe exci- 
tations of open strings attached to the p-brane. These open string vertices have the form 
£n 1 .../j, n dX >J ' 1 ■ ■ ■ dX^ n e lkX with certain Lorentz tensors £ Atl ... M „ and with momentum k par- 
allel to the Neumann directions. The case n = 1 (where the polarization e M is transversal) 
contains the massless modes: gauge fields living on the brane world-volume. 

The partition function of the theory with Neumann boundary conditions on one side and 
Dirichlet on the other follows from the boundary states as explained above, 

Z ND (q) = Tr HND q H ™ = {N\ q L °~& I D(x ) ) 



( 1 — r c V~^°° 1 — 

= —= (0|e~^n=i - a ™ a ™ q L o~24 eZ^ m =i | ) . (3) 

v2 

Orthonormality of the Fock ground states implies that only the contribution from the vac- 
uum sector survives in the second line. In particular, Znd is independent of the parameter 
xo. Computation of the vacuum expectation value above is straightforward. The result 
can be written as 

00 1 00 

ZND(q) = q^H(l-q n - l >)- 1 = -r- ) E^"^ • ( 4 ) 

n=l VVQ) n= i 



4 



Our main conclusion concerns the conformal weights of the boundary fields that can induce 
a transition between Dirichlet and Neumann type boundary conditions. The lowest weight 
that appears is h = j^. Above this value, the spectrum of conformal weights has half- 
integer spacings. The boundary condition changing operator with conformal weight h = jg 
corresponding to the lowest-energy state |cr) in the whole sector Hnd will be called a(x). 
We will also refer to a(x) as a "twist field" since the sum of irreducible Virasoro characters 
in (4) can alternatively be regarded as the character of a twisted U(l) representation. The 
absence of a vacuum state and the half-integer energy grading are symptoms for the fact 
that the jump from Neumann to Dirichlet destroys the simple U(l) Ward identities that 
are present in a boundary CFT with constant Neumann or Dirichlet condition all along 
the boundary. See [21] for general results about twist fields and partition functions in 
boundary conformal field theory. 

It will be our main concern in the following to find "substitutes" for the broken Ward 
identities, namely twisted Knizhnik-Zamolodchikov equations. 



2.2 Construction of the basic fields. In order to construct a field X(z, z) obeying the 
boundary conditions (1) on H = Hnd, we introduce a set of oscillator modes a r labelled by 
half-integers r G Z+ |. (We drop the superscript ^ for operators of the upper half-plane 
theory.) They are supposed to obey the relations 

[a r ,a s ] = r5 r - s , a* = a_ r . 

Creation operators a r , r < 0, generate the Fock space TC out of the ground state |cr), which 
is annihilated by modes a r with index r > 0. All the fields we shall consider act on this 
state space Ti. It is simple to verify that the decomposition X(z, z) = X{z) — X{z) of the 
bosonic field yields the desired properties if 

X{z) = x + % — z ~ r > = i E ~z~ r . 

To make the square root well defined, we have to introduce a branch cut in the plane, which 
extends from x = to — oo. Once the bosonic field is known, we obtain chiral currents as 

J(z) := idX(z,z)= ^ a r z~ r ~ x , J(z) := idX(z,z)= — ^ a r z~ r ~ x . 

Finally, the components T(z) and T(z) of the stress energy tensor are given by 

T(z) = lim - (j(w)J(z) - - — ~ ttt J , 

and likewise for T(z). Since T and T are quadratic in J and J, they satisfy the usual 
boundary condition T(z) = T(z) all along the real line Imz = 0. By the usual arguments 
[8] this implies that the modes 
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obey commutation relations of the Virasoro algebra with central charge c = 1. Here 
C_|_ — C- is a closed oriented contour surrounding the origin, with C+ contained in the 
upper half-plane and C_ contained in the lower half-plane. The commutation relation 
between L n and a r is easily checked to be of the form 

[ L n , Cl r ] V Cl n + r . 

It is convenient to introduce two generating fields T and J by the formal sums 



T(W) = L n ™- n ~ 2 , JH = E 



-1 — 1 

a r w 



One may think of T as being defined on the entire complex plane with T(w) = T(w) in the 
upper half-plane, Imw > 0, and T(w) = T(w), for all w with Imu; < 0. The generating field 
J naturally lives on the two- fold branched cover of the complex plane defined by u 2 = w. 
By introducing the branch cut from x = to — oo we have specified a coordinate patch on 
this surface with the local coordinate being denoted by w. In this chart, J(w) = J(w) for 
Imw > and J(w) = —J(w) for Imw < 0. 

The commutation relations for the modes L n ,a r with the bosonic field X(z,z) can be 
expressed in terms of T and J as follows 

[ T(w) , X(z, z) } = dX(z, z) 8{z - w) + dX(z, z) 8{z - w) , 
[ J(iu) , X(z, z) ] = % S(z — w) + iS(z-w) , 



where 
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We state two simple consequences of these formulas that shall be important below. We 
split J and T into two parts J(w) = J>(«0 + ^ < (w) and T(w) = T y (w) + T < (w) such that 



■l>(w) := a r w r 1 , T > (w) := L 

r>l/2 n>-l 



In the next subsection, we will use the commutation relations between the singular parts 
T>, J> of the generating fields T, J and the bosonic field X(z, z): 

[ J>H , X{z) ] = - (^) V2 -i- , [T>H , X(z)} = d z X(z) . (5) 

Moreover, we will need the following lemma, a proof of which is given in Appendix A. 

Lemma 1: One may rewrite the generating field T(w) in terms of the objects J>(«0 and 
J < (w), namely 

T(w) = 1(j < hjh + jhj > h) + -^^ . (6) 

z lb w 
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2.3 Bulk and boundary primary fields. Our next aim is to construct primary bulk 
and boundary fields. Here, the latter term refers to open string vertex operators which 
can be inserted on M < o or R>q without changing the boundary condition. They are in 
one-to-one correspondence to states in Hdd or Hnn, not to states in Hnd- We will see 
that bulk fields can be regarded as products of such boundary fields ^ g (z). Therefore we 
discuss these "chiral fields" first - admitting arbitrary complex insertion points, not just 
z E <9H. The fields \P g (z) are labeled by a real parameter g and enjoy the properties 

[T(w) , ^g(z)} = d z V g (z) S(x -w) + h tf g (z) d z S(z-w) , 

[JW»*g( z )l = g*g(*) S(z-w) . 

We have used the definition of the formal 5-function specified above and h — ^g 2 . For the 
commutators of the T > (w), J>(w) with the fields ^ g (x), this implies 

[T>(u;), * g (z)] = d z V g {z) + ^-A_ * g (z) , (7) 

[J>H,* g (,)] = (^) 1/2 ^^W • (8) 

Lemma 2: The unique solution (up to normalization) ^ g (z) to the requirements (7,8) 
is given by 

tf g ( Z ) = (i\ z ~h jgX^z) e igX > (z) where x> ( z j = ^ 



and X<(z) = X(z) — X > (z). Note that VP g (z) is normal- ordered, i.e., the annihilators 
a r ,r > 0, appear to the right of the creation operators. 

A proof can be found in Appendix A. 

Our next aim is to describe the U(l)-primary bulk fields (j)g(z, z). By definition, they obey 
the following commutation relations with respect to J> and T>, 

1 h 

i T >( w ) , <f>s(z, z) ] = d(j) g (z, z) + (f> s (z, z) (9) 

w — z [w — zy 

H — d(p g (z, z) + - — (f> s (z, z) , 

w — z (w — z) z 

[j>H, = 1/2 ^z- z - Q 1/2 ^M^~z) ■ (io) 

Note that each term from eqs. (7,8) appears a second time with z being replaced by the 
variable z. One can easily work out commutation relations between the full generating 
elements T(w), J(w) and the bulk primary fields <j> s (z, z). It is obvious from our discussion 
of boundary fields that bulk fields <f> s (z,z) can be written as products of chiral vertex 
operators, 

<f> g {z,z) = q g (z)q- g (z) . 



The formulas we have reviewed here would enable us to perform a direct computation of 
arbitrary correlations functions G(z) for bulk- fields <j) s (z, z) with two twist fields a inserted 
at x = and x = oo, 



G(z) := (a\ 0i(zi,zi) • • .(f) n (z n ,z n )\a) where <f> v (z v ,z v ) = ^ gv (z v , z v ) . 



The calculations would proceed by moving all annihilation operators to the right until they 
act on the ground state |<r). The same techniques would apply if there are extra boundary 
fields VPg(x) inserted in addition to the bulk fields <f> s (z, z). Since we will develop another, 
more generally applicable approach to the computation of correlation functions below, we 
do not enter details here. 



Before we conclude this subsection, we would like to derive bulk-boundary operator product 
expansions which allow to expand our bulk fields (j) g (z, z) in terms of boundary operators 
[10]. The essential idea is that, in the presence of a boundary, bulk fields split into products 
of chiral vertex operators inserted at points which are obtained from each other by reflection 
at the real axis and with opposite charges ("method of image charges"). In order to obtain 
concrete formulas, we first rewrite the bulk fields (j) g (z,z) = fy g (z)^ - g (z) in terms of 
X> (z,z) = X>(z) - X> (z) and X< (z,z)= X< (z) -X<(z), 

. x 2h / r- , /=\ 2h 

i\ , / \/z + Vz 



</> g (z,z)= ^-j &)~ h elgX<iZ,S) ^ X>{Z ~ Z) ■ (11) 

We have used the expression in Lemma 2 and then normal-ordered the right hand side 
with the help of the BCH formula, which leads to the additional y/z- and y/l- dependent 
factor. 



Lemma 3 (bulk-boundary OPE): For arguments z = x + iy close to the boundary, i.e., 
y > small, the operators <f>g(z, z) can be expanded in a series involving boundary primary 
fields, with leading asymptotics 



g*gX0 

0g( z > f ) ~ -^Ih 1 for x > , 

(j) e (z, z) ~ y 2h *2 g (^) for x<0 , 



where 1 is the identity field. 

Proof: Let us begin with the case x > in which y/z — — »• as y becomes very small. 

( i\ 2h 2h f2x\ 2h 1 



8 



We have also used that X K (x,x) = x and X > (x,x) = for x > 0, which is a direct 
consequence of the Dirichlet boundary condition. 

If x < and y tends to zero, the sum y/z + y/z" vanishes and we can estimate the behaviour 
of (f) g (z, z) according to 



(. x 2h / - \2h 

M {z - z) - h ( *"* ) e **<<*.*) e **><*.*) 
£J \Z + Z — Z\/ZZ/ 

(■ \2h / ■ \2h 



Observe that boundary condition changing operators themselves do not arise from the 
bulk-boundary OPE of bulk fields. Let us finally note that operator product expansions 
for bulk fields or for boundary fields can be worked out with the same techniques. Those 
for bulk fields <j> g (z, z) of course agree with the usual OPE of primary fields in the bulk. 

2.4 Correlation functions and the Knizhnik-Zamolodchikov equation. This 
subsection contains the main result of this section, namely a derivation of the Knizhnik- 
Zamolodchikov equation for correlation functions of bulk and boundary primaries in the 
presence of a transition from Dirichlet to Neumann boundary conditions at the origin. Let 
us look more closely at correlation functions containing n chiral fields VP g (z)> 

F(z) := (a\^ 1 (z 1 )...^ n (z n )\a) where V v (z v ) := * gl/ (z„) • 

Before we start, let us state two elementary formulas for the action of J>(w), T > (w) on 
the ground state \a) : 

T>H \a) = (^h a + -L_A \a) , (12) 
J > (« ; ) \ a ) = a r \a)w- r - 1 = , (13) 

r>l/2 

where h a is the conformal weight of the state |a). We will recover the equation h a = 
in a moment. The object (cr| dual to \a) obeys the relations (a\ J<(tu) = = {a\ T < (w). 

A first differential equation is obtained by inserting the generating field T(w) into the 
correlation function: 

(a\T(w)^ 1 (z 1 )...y n (z n )\a} = (a\ J > (w)^ 1 (z 1 ) . . . ^ n (z n ) \a) 



< * \ 11) 7 



V = \ 



w — z v " ' (w — z v ) 2 ) w 2 



F(z) + - (a\y 1 (z l )...y n (z n )L_ l \a) . 



w 



Here, we have commuted T>(w) through the fields ^ v (z v ) until it acts on the ground state 
| a) so that we can use formula (12). 
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Now we want to compute the same correlation function with the help of the affine Sugawara 
construction, i.e. by exploiting eq. (6): 



(a\T(w)q 1 (z 1 )...q n (z n ) \a) 



1 1 

= - (a\ J > (w)J > (w)^> 1 (z 1 ) . . . ^n(Zn) W) + <crj . . . V n (z n ) \a) 



v=l 



1/2 H J>H*i(zi) • • .*„(*„) W) + f(z) 

w — z v 16w z 



\JZi/Z^ 



+ 



w (w — z v ){w — Zfj,) 16w 2 



F(z) 



Comparison with our first formula for the insertion of T(w) yields h a = j^. From the 
residue at w = we get 

(<r\9 1 (z 1 )...9 n (z n )L- 1 \v) = \ J2^= F ^ • 
Finally, from the residue at w = z v we obtain 



d Zv F{z) = 



hiL. _|_ ^2 I— 

Zv , V Zv Z v Zn 



F(z) . 



(14) 



This is the Knizhnik-Zamolodchikov equation we were after. Note that the terms in square 
brackets determine a flat connection, as in the ordinary Knizhnik-Zamolodchikov equation. 

We can solve (14) by a simple coordinate transformation. In fact, if we introduce coordi- 
nates Uv = \fz v and the function F u {u\, . . . , u n ) = Y[ u u^F(uf 7 . . . , u^), then the system 
(14) of first order differential equations becomes 

(2 n \ 



This equation is formally identical to the usual Knizhnik-Zamolodchikov equation with 2n 
fields of charges ±g u inserted at the points ±u v . The solution to (15) is given by 



v=l 



l<v<fj,<n 



Uv + U^- 



(16) 



The free parameter k can be determined from the boundary condition of the bosonic field 
X on the positive real line, i.e. k = «(xq) depends on the position xq of the D-brane. 
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3. Twisted Knizhnik-Zamolodchikov equation for multiple 
transitions 



In the following, we study n-point functions of a free bosonic field theory on the half-plane 
with several insertions of twist operators placed along the boundary. In the corresponding 
string diagrams, open strings stretch between three or more branes of various dimensions. 
As long as only one DN-jump occurs, a simple Hilbert space formulation of the boundary 
CFT is available, and we can solve, in principle, for correlation functions by purely algebraic 
techniques, as indicated in the last section. In the presence of many boundary condition 
changing twist fields, it may be simpler to resort to OPE methods and to the theory of 
complex functions on higher genus Riemann surfaces, and this is the approach we pursue 
in the present section. We begin with a very brief review of relevant input from the 
theory of hyperelliptic surfaces. We then discuss Ward identities in the second subsection. 
The latter allow us to derive a system of linear first order differential equations for the 
correlation functions similar to the Knizhnik-Zamolodchikov equations. Note that free 
bosons on higher genus surfaces without boundaries (i.e. higher loop diagrams of closed 
strings propagating in a flat target) have been studied in great detail in [4,31]. 

3.1 Hyperelliptic surfaces. Our aim is to investigate a scenario in which a bosonic 
field X(z, z) is defined on the upper half-plane with boundary conditions switching between 
Dirichlet and Neumann at 2g + 2 points Xi, i = 1, . . . , 2g + 2, on the boundary. Without 
loss of generality, we shall assume that X2 g +2 = oo =: xq- To be more precise, we impose 
Dirichlet boundary conditions in the intervals ]xi, for % odd and Neumann boundary 
conditions along the rest of the boundary, i.e., 

X(z,z) = Xq for z = z e D k :=]x2k-i,X2k[ 

and 

dyX(z,z) =0 for z = z E N k := }x 2 k-2, x 2k -i[ ■ 

The variable y is defined through z = x + iy, and fc=l,...,<7+l. In terms of the chiral 
currents J{z) = idX(z,z) and J(z) = idX(z,z), these conditions become 

J(x) = —J(x) for x G Dk , and J{x) = J(x) for x £ N k 

As in the previous section, it is convenient to work with a single field J that contains all 
information about the two chiral currents J and J. Such a field necessarily lives on a 
two-fold branched cover of the complex w-plane, namely on the hyperelliptic surface, M, 
of genus g which is described by the equation 

2g+l 

uj 2 = P( w ) : = Y[( w - Xi ) . 

i=i 

Introducing branch cuts along the intervals N k = \x 2k -2i #2fc-i] > we obtain a particular 
coordinate patch of this surface with local coordinate w. In this chart, -i(w) satisfies 
J (to) = J(w) for Imw > and J(w) = —J(w) for Imw < 0. The Virasoro field T obeys the 
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gluing condition T(x) = T(x), all along the boundary, since it is quadratic in the currents. 
Consequently, the generating field T(w) is defined on the complex tu-plane and coincides 
with T (resp. T) on the upper (resp. lower) half-plane. 

The coordinate w on the complex plane lifts to a meromorphic function of degree 2, also 
denoted by w, on the hyperelliptic surface M. This function defines a two-fold covering of 
the sphere branched over 2g + 2 points Qi, . . . , Q2 g +2, where w(Qi) = X{. A basis for the 
space of holomorphic 1-forms on M is then given by 

w k ~ 1 dw . , 
V p ( w ) 

It will be convenient to work with a canonical homology basis {7^, %} on M chosen as in 
Figure 2. We denote by Qki the period of the 1-form ui along the cycle 7^, i.e., 

/ / w l ~ 1 dw 



7fe J 1k 



The basis of holomorphic 1-forms, {Cfc}, dual to the canonical homology basis {7/c 7/0} is 
defined by the equation 

Ci = Ski ■ 



In particular, we have the relation 



7fc 



1=1 



and the matrix O is invertible. The period matrix r is given by 

Tkl ■= j> Cl , 

•J Ik 

and it is known to be symmetric and to have positive definite imaginary part. The surface 
M has an anti- holomorphic involution induced by complex conjugation on the complex 
plane. In terms of the functions w and 00 = ^P{w), it can be written as 

(w, ui) — > (tD, u) . 

This involution will be used to extend the theory from the upper half-plane to the lower 
half-plane while taking care of the boundary conditions on the real axis. There is a second 
holomorphic involution that interchanges the two sheets of M and that can be written as 

(to, u) — > (w, —uj) . 

This involution is used when passing from the sphere with cuts to its cover M. 
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3.2 The Ward identities. To begin with, we introduce the correlations that we 
plan to investigate below. Besides the primary bulk fields (j)g(z, z) := exp(igX (z, z)), they 
involve 2g+2 boundary twist fields inserted at the points Xi, which induce changes between 
Dirichlet and Neumann boundary conditions. From our discussion in the previous section 
we know that there is an infinite number of boundary condition changing operators that 
we could insert. But it is sufficient to study the fields a(x) of conformal weight h = jq, 
since all others can be obtained out of a(x) by OPE with chiral fields. Thus, our discussion 
deals with correlators of the form 

G(z,x) = (0i (zi, zi) . . .(j) n (z n ,z n ) cr(xi) . . .<j(x 2g +i)) (17) 

where we use the notation <p v = 4>% v i an d where the boundary field (j(x2 g +2) = cr(oo) is 
absorbed in the notation (...) = (a\ . . . |0), with |0) denoting the vacuum state. We could 
insert further boundary fields VP g (z). Their interpretation depends on whether they are 
inserted in an interval with Dirichlet or Neumann boundary conditions. In the former 
case, they could induce jumps in the Dirichlet parameters Xq if they are associated with 
open strings stretching between branes at different positions. A boundary operator ^ g (z) 
inserted in one of the Neumann intervals, on the other hand, creates an open string which 
has both ends on the same brane (a Euclidean 1-brane, in our case) and moves with some 
definite momentum along its world-volume. 

As far as Ward identities are concerned, correlation functions of such boundary fields are 
actually more fundamental, since one may split the bulk fields <f> s (z, z) into a product of 
\Pg(z) and ty-g(z). For this reason, most of our investigations below involve the correlation 
functions 

F(z,x) = (tfi(zi)... a(x 1 )...a(x 2g+ i)) , (18) 
from which the correlators G(z,x) can be reconstructed. 

Our analysis will make essential use of the Mittag-Leffler theorem, and hence it is based on 
the study of singularities in correlation functions. The latter are encoded in the operator 
product expansions between chiral fields and the bulk and boundary fields appearing in 
(17,18). For the Virasoro field T one has the standard expansions: 



T(w)* s (z) 
T(w) a(x) 



h 



(w — z) 2 
(w — x) 2 



+ 



1 



d z 



+ 



w — z 

1 



w 



X 



d x 



*g(*) 
a{x) , 



(19) 
(20) 



Here and in the following, the symbol ~ means "equal up to terms which are regular as 
w — > z" . According to the rule <f> s (z,z) ~ ^ g (z)^- g (z), the operator product expansions 
of T with (j) g contain further terms in which z is replaced by z (note that h g = h- g ). These 
formulas may be compared with equations (7) and (12) in Section 2. For our correlation 
functions, eqs. (19,20) imply 

(T(w) i(zi) . . . y n (zn) <r(xi) . . . (T(x 2 g+i)) 



d 



~\(w- z v ) 2 w - z v dz v 



2g+l 

E 

i=l 



d 



(w — Xi) 2 w — Xi dxi 



F{z,x) .(21) 
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The situation is more subtle for the current J(w), which we recall is only well-defined on 
a surface of genus g. More precisely, J(w)dw is a meromorphic 1-form on a hyperelliptic 
surface. For the operator product expansion of -l(w) with the field ^ s (z), we shall use 



J(w) tyJz) dw 



w 



VJz) dw 



(22) 



Indeed, the right hand side has a first order pole at w = z with residue g and is regular 
otherwise. Equation (22) generalizes formula (8) in Section 2.2. To determine the oper- 
ator product expansion between J and the twist field cr(x), we observe that the leading 
contribution 

J(w)a(x) ~ (w - z)^" 1 "^ t(x) + . . . 

must involve a field r of conformal weight h T = h a + 1/2 + Z. Otherwise, the expansion 
would not be consistent with the periodicity properties of J close to the branch point 
at w = x. Among the boundary condition changing operators, there is one field with 
conformal weight h T = \ + jq which gives rise to the most singular contribution diverging 
with (w — x) -1 / 2 , cf. the spectrum (4) computed above. 

After multiplication of the previous equation with dw, we are supposed to study the right 
hand side in the local coordinate £ = \/w — x. The outcome is rather simple: The form 
(w — x)~ 1 / 2 dw = 2d£ on the right hand side is regular at £ = so that we conclude 



](w)a(x)dw ~ 



(23) 



i.e., the singular part of the operator product expansion between ]{w)dw and the twist 
field a(x) vanishes. 

As we will see, the following important formula is a consequence of these operator product 
expansions: 

(JH . . . ^ n {z n ) a{ Xl ) . . . a(x 2g+1 )) 



where 



E 



ak(z,x)w k 1 



F(z,x) , 



(z,x) = ^n^ 1 zA fc x -^g,S fc (^) 



1=1 



v=l 



(24) 



and the parameters A fc x := Xq — Xq +1 are obtained from the values of the bosonic field 
at the boundary. 

To derive the formula (24) we exploit the fact that a meromorphic 1-form on a compact 
Riemann surface is determined by its principal part up to some holomorphic 1-from. Op- 
erator product expansions, on the other hand, contain all information about the principal 
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part. Hence, from eqs. (22,23) we conclude that 



(JH i(zi) . . . V n (z n ) a( Xl ) . . . a(x 2 g+i)) 

Note that the insertion points z v and parametrize a whole family of meromorphic 1- 
forms, and the coefficients (3 k may depend on them. Actually, we can determine this 
dependence completely. To this end we integrate the above equation along a loop 7^ which 
surrounds the interval [x2k, #2fc+i]> as shown in Figure 2. On the right hand side of our 
equation, this integral may be expressed in terms of the matrix O and 



B\ Z ) := / -=L= dw . (26) 

Note that the matrix elements Q k i and the functions B k (z) depend on the insertion points 
Xi. With these conventions we find 

P n g 

f dw (r.h.s. of (25)) = Sfc (^) F (*> f ) + E ^ A • 

Next, let us analyze the integral over the left hand side of equation (25). By a deformation, 
we can make the integration contour 7^ symmetric under a reflection 7^7 along the 
real line. Now, we may split 7^ into two parts 7^ , 7^ with the property Inry^ > and 
Im7^ < 0, so that each piece lies entirely in one of the half-planes. With these conventions, 
our contour can be written as a composition 7^ = 7^ o 7^ which obeys 7^ = ~l k • ^ we 
recall, in addition, that the field J coincides with J on the upper and with — J on the lower 
half-plane we deduce 

j> J(w) dw = / J(w)dw+ / J(w)dw = i / dX(w,w) = z(xq-Xq +1 ) . 

In the penultimate step, we have expressed the currents through the bosonic field X by 
J(w) = idX(w,w) and J(w) = i dX(w,w). The contour integral over the differential dX 
is finally determined by the values of X at the boundary. For the integration over the left 
hand side of (25), this result implies that 



I dw (l.h.s. of (25)) = i(x*-^ +1 )F(z,x) =: iA k x F(z,x) . 
J ik 



Putting all this together, we arrive at the following formula for the function (3 k (z,x): 
p k (z,x) = J2 Q kl' (iA'xo-^g^ 1 ^)) F(z,x) =: a k (z,x) F(z,x) . 



1=1 \ v=l 
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The functions introduced here depend on the insertion points z„,Xi, the charges g u and 
the values Xq of the bosonic field at the boundary. Additional information, e.g., on the 
unknown function F(z,x), is not needed. This concludes our derivation of eq. (24). 




Figure 2: The curves 7^ run counterclockwise around the Neumann cuts on the 
real line. The curves 7^ run clockwise and close on the second sheet through 
Neumann intervals. 



3.3 The twisted Knizhnik-Zamolodchikov equations. We now start to exploit 
formula (24) together with the Sugawara construction of the energy-momentum tensor T 
to compute the effect of inserting the field T into our correlation functions. If we include 
one extra field ^ g (u) into eq. (24), differentiate with respect to it and set g = 0, we obtain 
as a consequence of J (it) = -d u ^ g (u) | g= o 

(JH J(u) l(zi) . . . ^n(Zn) 0{xi) . . . <T(x 2ff +l)> 



y - 



9 n 



it - z,, 



a k {z,x)w k 1 g u y/P(z u ) 



=^1 v^Rv^H 



10 - z. 



+ E 



;; akaiw^u 1 - 1 , 1 d / VP(u) 



+ 



1 D^„.\„..fe-1 



^ y/P(w)y/P(u) /PH * 



lt 



F(z, x) 



At this stage we can subtract the term 1/(10 — w) 2 , multiply by a factor 1/2 and perform 
the limit u — > k;. A short and elementary computation gives 



(T(k;) *i(zi) . . . ^ n (^ n ) (t(xi) . . . <t(x 2 «h-i)) 
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1 gi/g/x V P M V P M { s^^ ak(z,x)w k 1 g u y/P(z u ) 



2 ^ P(w) w - z v w - 



fc=l V=\ 



P(w) 



w - z. 



1 



+ E 



a k uiw 



k+l ~ 2 i q^p'H'^" 1 



P(tu) 



4y / PH 2 



E 



k,l=l v ' " V ' V"V k,l=l 

The same correlator has been computed in eq. (21) directly with the help of operator 
product expansions between the Virasoro field T and \I/g, a from (19,20). Comparison 
of the residues at w = z v in the two different expressions gives the ^-components of the 
Knizhnik-Zamolodchikov equations: 



d Zv F(z,x) 



E 



— 7... — T, ^ 



gvOL k z; 



fc-1 



F(z,£) 



If we restrict to the case of two twist field insertions, i.e. (7 = 0, the last term vanishes. 
Putting xi = and hence P{w) = w, formula (14) is recovered. 

Computation of the residues at w = Xi yields a formula for the derivative of F with respect 
to the position Xi of twist fields. Using that 



Res 



k,i=i 




^ l 1 B l (w)'w k - 1 



« Z-^i r f .._r r . 



E <v 

fc,Z=l 



71 



1 X?" 1 J_ 

2 v^O^-x, 



Egv^^) , /-> ->\ fc-i 



- 1 - T^— -- V o- 1 

8 ^ Xi - Xi 4 ^ kl 

jjti 3 k,l=l 



k=l 

x k ~ l 
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F{z,x) .(27) 



To conclude this section, we come back to the original correlators G(z, x) of bulk-fields 
(f) g (z,z) and boundary twist fields a(xi). The differential equations they obey will be 
formulated with the help of the following functions: 

>W \fp~W) 



u (z,w) = 
Xo(z) = 
o-a(z) = 



1 



z — w 



1 23+1 I 

_1 y _J__ 

9 £-~> 7 — T: 



z — w 



9 Jfc-l 

V - 



0- 1 (^H-S'(«;)) , 



1 ^/Pjt) A 



i=l 



Z Xi 



^-E 



z — z 



k,l 



n;; 1 (£'(*) -B'(z)) , 



y -4— ^ A'x 
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Taking into account the equation (f> g (z, z) = ^ g (z) *&- g (z), we conclude that the correlation 
functions G(z, x) must satisfy the following set of first order linear differential equations: 

d^G(z,x) = A^G(z,x) for all £ = z u ,z u ,Xi (28) 

with the connection matrices A Zv , A Zv and A Xi being defined by 

A Zv = ^2gvgpWo(z v ,Zp) + gt,\o(z v ) + ig u <TA(z v ) , (29) 
A Sv = J^gi/g^ota/,^) + gl\ (z„) +ig„(TA(zu) , (30) 



i f n V 

A Xi = - lim (x-Xi) [y^g v uo(x, z v ) + ia A (x)\ 
2 ™ \ti J 



8H *M 4 k% J-n y/PGKt-*i) 

We have introduced the function Hi(xi) := Ylj-^^Xi —Xj). Note that the term in brackets 
has a simple pole at x = Xi which we cancel by the extra factor x — Xi before performing 
the limit. Equations (28) through (31) constitute the main result of this section. 



4. Construction of correlation functions 

It remains to reconstruct the correlation functions G(z, x) from the system of linear first 
order differential equations that we obtained in the previous section. Integration of the 
equations is, in principle, straightforward, but it leaves one constant factor undetermined. 
The latter is found explicitly in terms of the boundary conditions. Moreover, we shall 
manage to express the correlators G(z, x) in terms of rather elementary building blocks. 



4.1 Integration of the z-connection. To begin with, we simplify our problem by 
fixing the insertion points of the boundary twist fields a(xi) and considering only the 
dependence of Gg{z) = G(z,x) on the positions (z u ,z v ) of bulk fields. This means that 
we have to integrate the ^-connection A Zv dz v + A Zv dz v defined in eqs. (29,30) . The result 
will be written with the help of two functions Gq(z, w) and Sq(w) which are given by 



G (z,w) := 2 Re / cu (£, w) 

J X 



S (z) 



lim 

V— >x 



2 Re j A (f ) d£ - \og(v - v) - Re log P(v) 

J V 



(32) 
(33) 



Here, the point x is chosen to lie in the Dirichlet-interval ]xi,X2[. The integrand in eq. 
(32) is regular on the real axis and hence the integral is well defined. One can easily see 
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that its value neither depends on the starting point x G ]xi,X2[ nor on the choice of the 
curve 7 from x to z. In contrast, the integrand Ao(£) in our definition of Sq(z) has a simple 
pole on the real axis. This is the reason why we subtract the divergent term log(i> — v) 
before taking the limit v — * x. The contribution — Relog-P(i>) is added to render the whole 
function independent of the integration contour and, in particular, of the starting point x. 
More details on the definition of Sq(z) and a discussion of its properties can be found in 
Appendix B. 

The two functions we have just defined possess a number of abstract properties that char- 
acterize them uniquely. First of all, it is not difficult to see that Gq(z,w) is simply a 
Green's function on the upper half-plane, i.e., it obeys 

A z Gq(z, w) = 4nS(z — w) for Imz > 

and the boundary conditions 

d 

G (z,w) = for z G Di , G (z,w) =0 for z G N t . 

olmz 

A standard computation shows that Go(z,w) is symmetric in its arguments, Gq(z, w) = 
G (w,z). 

The function So (z) , on the other hand, is harmonic throughout the whole upper half-plane, 
i.e. A z So(z) = 0. It diverges at the boundary with a leading singularity of the form 

S (z) ~ =F log \z — z\ + . . . for Re z G <j ^ . (34) 

We are now in a position to integrate the differential equations (29,30) for the correlators 
of bulk fields <j> s (z, z). The result is given by 

log^(z) = [jA Zv d^ + A- z J^) + Hw) 

J w 

1 n n n g+1 „ ^ 

= g E g^GoC^ + X^SokO + E^ E / fc^^' ( 35 ) 

",M=1 !/=l I/=l i = l ^ 



In the first line we have chosen some arbitrary curve in the configuration space of n particles 
in the upper half-plane starting at points w v with ImtOj, > 0. The integration "constant" 
A(w) depends on the choice of the starting point and has to be fixed such that the resulting 
function \ogG$(z) satisfies the desired boundary conditions. In passing to the second line, 
we have extended the integration to w v = x on the boundary and inserted the definitions 
(32,33) of the functions Go and So- Then we use the auxiliary formula 
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To prove this formula one should notice that the function on the l.h.s. of the equation is 
harmonic in the upper half-plane, that it satisfies Neumann boundary conditions along 
the intervals Nk and that it approaches the constant values Xq for z G Dj~. By explicit 
computation, one can establish the same behaviour for the function on the r.h.s. Since 
these properties are sufficient to determine the functions uniquely, the desired equation 
follows. We employ it to bring the third term of eq. (35) into a form that allows the most 
explicit control of the boundary behaviour and hence is quite appropriate for fixing the 
remaining A (a;). 



4.2 Integration of the x-connection. Before we address the integration of the full 
connection (29,30,31) below, we investigate another simplified situation in which there 
are no bulk fields present. Consequently, the charges g v in the equation (27) can be set 
to zero, and we are confronted with the problem of solving the following equation for 

Z A (x) :=G(x): 

Vfc=i / 4 2 m=i ^ v p (0(£-O 

where Hi denotes the function Hi = Ylj^i( x i ~ x j) as before, and A = {A^x }. These 
differential equations were solved by Zamolodchikov in [42]. Here we simply quote the final 
result: 

2g+l 

Z A (x) = n^-^O'^det-^^e^^M-^ ^" ^ . 

i>j 



4.3 The correlation functions G(z,x). The results of the previous two subsections 
can be combined into explicit expressions for the correlators G(z, x) of bulk fields 4> g (z, z) 
and boundary twist fields a{x). To see this we note that solutions of the differential 
equations (28) can be found by integrating the connection 1-form A Wv dw u + Ayj v dw v + 
A^d^i along an arbitrary curve 7(f) = (7^ (t), 7x(t)), t G [0, 1], that ends at the point (z, x) 
in the (2n + 2g + l)-dimensional real configuration space. 

With some care (see the first subsection) we can start the integration with all insertion 
points being on the real axis. Furthermore, we may choose 7 such that all twist-fields are 
moved to their final position at x before we begin moving the bulk fields into their desired 
locations. 

In more mathematical terms this means that 7 consists of two parts 7 = 7^) o 7W with 
da { }\t G [0, §]) =0 and d t ^f\t G [§, 1]) = 0. As long as t < J, the bulk fields are 
located at points w v = 7^(t) belonging to the first Dirichlet interval D\ =]xi,X2[. This 
implies u>o(x, 7^(t)) = for t G [0, |], so that one term in our expression (31) for A x . 
drops out. Hence, we are precisely in the situation considered in the previous subsection, 
and the integration of our connection 1-form over 7 in the interval t G [0, -|] gives Z A (x). 
When we continue the integration to t = 1, we add the expression for logG^(^) computed 
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in eq. (35). After some rewriting, the following result for the correlations function G(z, x) 
is obtained: 

G(z,x) = Z A (x) exp^g^ (^)j exp^zg^S- 1 ^^)^ • ( 36 ) 

Here, ^ {zv) denotes the potential that is created by n — 1 charges g M at points z M 7^ z v 
in the presence of the boundary with mixed boundary conditions, 

1 f 1 9+1 f X2i 8 

= / d 2 wG Q (z l ^w)J2s,S(w-z,) + - / ^— - -Go&^dt . 

It is quite instructive to interpret each of the three factors in our final expression for the 
correlation function G(z, x) directly within conformal field theory. For the moment, let 
us specify the number of bulk fields in G by some extra superscript, i.e. we shall write 
G{z, x) = G^ n \z, x). Now consider the object 

~, , 1 , ( G^(z,x) 
${z v ) := —log ' 



ig„ ta V exp(g 2 ^o(^))G'(- 1 )(f',x) 

where z' denotes the set of n — 1 bulk coordinates z^, \i 7^ v. It is easy to determine the 
behaviour of $(z u ) as a function of the bulk coordinate z v from the bulk and bulk-boundary 
operator product expansions of the fields (j)g{z v , z v ), cf. Subsect. 2.3: 

A Zv $(z v ) = 2n 6 ( z » ~ z ») 

$(z v ) = xj for z v e Di , — $(z v ) =0 for z e . 

olmzv 

These properties characterize the function $(z v ) uniquely, and by standard formulas from 
electrostatics we obtain that &(z u ) = ^^^(z^). An iteration of this construction along 
with G(°\x) = Z A (x) leads to our product formula (36) for the correlation function 

G(z, x). 

The three factors can be interpreted as follows: Z A {x) is a "partition function" corre- 
sponding to some line charge distribution provided by the twist fields alone; the term 

n 

is the electrostatic potential corresponding to a configuration of n point particles with 
charges g\, . . . ,g n located at the points zi, . . . , z n and line charges distributions along the 
Dirichlet intervals. The term 

n 

E dl S {z„) 
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can be interpreted as a renormalized electrostatic self-energy of the point particles located 
at Z\ , ... , z n . 

4.4 Path integral approach and extensions. Before we conclude, let us briefly 
sketch how the theory can be formulated in the path integral approach. This is important 
for the following two reasons: First, as long as one is only interested in free bosons, the 
path integral approach is a powerful alternative to our analysis above involving Knizhnik- 
Zamolodchikov connections. The path integral formulation allows for a more direct com- 
putation of the correlation functions (but it does not easily extend to non-abelian group 
targets.) Secondly, using path integrals we will be able to describe rather easily possible 
extensions of our analysis to compact targets and to D-branes with B-fields. 
To begin with, we consider once more the familiar situation of a non-compact 1-dimensional 
target and Dirichlet parameters Xq. We denote by G the Green's function of the Laplacian 
on the upper half-plane, 

A z G(z,w) = S(z — w) , 
subject to the boundary conditions 

G(x, w) = for x G D , 
d y G(x, w) = for x E N . 

The Gaussian measure with covariance G and mean is denoted by and we use x 
for the corresponding random variable. With the help of the Dirichlet parameters Xq, we 
define a real-valued function £ on the upper half-plane by 

A£ = 0, 

fll* = X » 9y^\ N . =0, I = . 

The effect of the Dirichlet parameters is incorporated through a shift of the random variable 
X by £ which gives us the bosonic field X = x + £• It appears when we construct the basic 
fields of the theory, namely the vertex operators 

(p g (z,z) = :e igX ^~ z) : . 

In this framework we could recover the correlation functions above by integrating products 
of fields (p g (z, z) using the Gaussian measure. 

When the free boson X takes values in a circle S 1 , the boundary conditions depend both 
on Dirichlet parameters Xq and on Neumann parameters denoted by y l . The Neumann 
parameters determine the Dirichlet parameters of the T-dual theory and can be thought 
of as the strength of constant Wilson lines turned on along a Neumann direction. In Sect. 
3, we have worked in a local chart with the Neumann intervals cut out. Information on 
Neumann parameters is, therefore, lost unless we take a second chart into account which 
has cuts along the Dirichlet intervals. As above, we can derive Knizhnik-Zamolodchikov 
equations for each chart; the full correlation functions of the compactified theory are to 
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be built up from the respective solutions in such a way that the boundary conditions are 
met. 

The path integral computation of correlators is rather easy to adjust to the compacti- 
fied situation: if we set Neumann parameters to zero and restrict attention to the fields 
<p g (z, z), we can use precisely the same formulas as above. 

In order to incorporate non- vanishing Neumann parameters yQ and to compute more gen- 
eral correlators for fields (p g c ^(z,z) with (g,g) taken from an even, self-dual Lorentzian 
lattice, we introduce a real-valued function rj on the upper half-plane defined by 

Ai] = , 

Vim = f , d y r]\ Di = , i = 1, . . . , g + 1 . 

Then, we set 

™ 9 A W ) := gG'(w,z) + i](w) , 

where G' denotes the Green's function of the Laplacian on the upper half-plane with 
interchanged Dirichlet and Neumann boundary conditions, i.e., 

AG'(w,z) = 6(w-z) , 
d y G'(x,z) =0 for x E D , 
G'(x,z) =0 for x G N . 

We now introduce the disorder operators D g (z,z) satisfying 

D g (z,z)F[dX] = F[dX + *dw g , z ] , 

for any functional F. The left- resp. right-moving chiral vertex operators can then be 
written as 

^g(z) = (pz(z,z)Dz(z,z) , 4>g{z) = ipi(z,z)D_s(z,z) ; 

see also [17] for more details and for an application to soliton quantization in 2-dimensional 
theories. The basic fields of the compactified theory are products 

where (g, g) lies in some even, self-dual Lorentzian lattice. For another approach to the 
rational compactified boson, the reader is referred to [22]. 

Another extension would involve the appearance of S-fields on our D-branes. This has 
attracted some interest recently, because of its relation with non-commutative geometry, 
see e.g. [15,35,36, 3,16] and references therein. Non-trivial B-fields can only exist if one of 
our branes is at least 2-dimensional. For simplicity, we shall focus on a pair of a Dp-and a 
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DO-brane. The field strength on the Dp-brane will be denoted by B. In terms of boundary 
conditions for a multi-component free bosonic field, the situation is described as follows 

d t X a (t,0) = and d a X a (t,n) = B a b d t X b (t,n) for a, 6=1, ...,p . 

The spectrum of the associated boundary condition changing operators and the Green's 
functions in the presence of two twist fields have been discussed at various places (see 
e.g. [36,11]). Our techniques from Sections 3 and 4 allow to extend such investigations to 
the case of multiple twist insertions. Instead of giving the details here, we simply state 
how one has to adjust the path integral computation to the new scenario. This is rather 
easy: All it requires is to replace the function G above by some matrix valued Green's 
function Gb = (G B b ). The latter is a Green's function for the Laplacian l p A on the 
upper half-plane (l p denotes the p-dimensional identity matrix), subject to the boundary 
conditions 

G a B b (x,w) =0 for x G D , 

d y G a B b (x, w) = iB a c d x G c i(x, w) for x G N . 
With the help of this function, the calculation of correlators proceeds as before. 

5. Outlook 

We have succeeded in decomposing the complete bulk and boundary correlators in the 
presence of DN-transitions into functions with rather natural interpretations - both from 
the point of view of electrostatics and from the CFT perspective. This is useful for carrying 
out the remaining step in the computation of string amplitudes, namely the integration 
over insertion points of fields on the world-sheet. The calculation of such string amplitudes 
gives effective actions involving a hyper-multiplet x which comes with the twist fields. To 
leading order, the bosonic part of these actions can be found in [14, 27,25]. Multiple twist 
insertions allow to compute higher order corrections. 

When we turn on a S-field, the string amplitudes may be described through field theories 
on some non-commutative space. It was suggested in [36] that these theories are related to 
some model on an ordinary commutative space through a complicated non-linear transfor- 
mation. This statement can be checked order by order in the effective description. After 
the appropriate (but straightforward) extension to non-vanishing S-fields, the considera- 
tions presented above may be used to perform a similar analysis for theories which contain 
a hyper-multiplet x- 

Keeping the bulk insertions fixed, the sequence of correlators with arbitrarily many twist 
field insertions can be viewed as building blocks of the perturbation series of a relevant 
perturbation by the twist field. This "tachyon condensation" is responsible, e.g., for the 
formation of D0-D2 bound states, as discussed in [23]. Upon integrating over twist field 
insertion points in the one-point functions Z&(x) exp{g 2 So(z)}, one would arrive at one- 
point functions which characterize the boundary theory after tachyon condensation. Sen's 
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approach [37] and the results of [33] allow one to circumvent the relevant boundary flow 
and to replace it by a combination of marginal bulk and boundary deformations. However, 
some questions as to the equivalence of both procedures remain open, and it might be useful 
to have an independent check of these methods. The correlation functions constructed here 
provide a starting point. 

For applications to superstring theory, it is mandatory to extend our analysis to free 
fermions. This does not pose serious problems, since systems of an even number of fermions 
can be bosonized. 

Problems of the type of our free boson problem are encountered in general boundary CFT 
as soon as the "parent" CFT on the plane admits different boundary conditions. For some 
general results on the rational case, see [20,21,22]. The spectrum of boundary condition 
changing operators can be derived as in Sect. 2, once boundary states for the "constant" 
boundary conditions are known. Again, the computation of correlators becomes non-trivial 
if boundary conditions with different gluing automorphisms are combined. In non-abelian 
WZW models, which constitute and important generalization of the free boson case, the 
Sugawara construction can be exploited in a similar fashion as for the free boson and 
leads to twisted, non-abelian Knizhnik-Zamolodchikov equations. The partition functions 
counting BCCOs in non-abelian boundary WZW models are linear combinations of the 
twining characters investigated in [18] (see also [5] and references therein). Apart from the 
models with affine Lie algebra symmetry, there is the rather large class of so-called "quasi- 
rational CFTs" [28] on the plane for which generalizations of Knizhnik-Zamolodchikov 
equations exist even without a Sugawara form for the energy-momentum tensor [2]. It 
might be interesting to see how such structures extend to boundary CFT. 
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Appendix 

A: Proofs of Lemma 1 and Lemma 2. 

Proof of Lemma 1 : We start from the usual expression for T in terms of J and rewrite it 
until we can perform the limit w± — > W2- 

JOi) J(w 2 ) - t 
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= J<(u;i)J(u;2) + J(tw 2 )J>(tWi) + [J>(«>i) , J(u>2)] 



(wi - w 2 ) 2 

i f^-Y /2 + \ (V) 1/2 - 1 
= J<(k;i)J(k;2) + J(w2)J>(wi) + Z 



We can now perform the limit w\ ^ W2 ='■ w to recover the generating field T(w) from 
the last formula 

TH = \{ J<HJH + JHJ>H ) + 

u l/2 +u -l/2_2 1 

where we used lim 



u-i 2(1 -w) 2 



Proof of Lemma 2: The derivation of the commutation relation with J>(w) is straightfor- 
ward. So, let us turn directly to the calculation of the commutator with T > (w). Recall 
that h = g 2 /2. Then, 

[T>(«;),* g (z)] 

= (j^J ([T>(«;),e lgX < (2) ]e lgX > (z) +e lgX < (2) [T>(«;),e JgX > (z) ]) 

e lsX < {z) (^ig[T > (w),X < (z)} - h [X K (x), [X < (x),T > (w)}\ 
+ ig [T> («;),*>(*)] + h [*>(*), [*>(*), T>H] ]) e lgX > « 



2z 



r,s>0,n> — 1 



r,s<0,n> — 1 
igX>(z) 



e 



1 


w 




z 




1 




w 




z 




1 




w 




z 



z 1 1 , ... , 

— 



\2z J V / — z) 2 wz 

d z VJz) H — + /i ( - — * ] # ( z ) 

8W w-z^ gw \{w-z) 2 z(w-z)J gv 7 

^* g (z) + ^A_*«( z ) • 

In the process of this computation we have inserted the commutation relation between 
L n , a r and eq. (5). The rest involves only standard algebraic manipulations. ■ 
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B: The function S Q (z). 

In this appendix we want to explain a number of properties of the function So(z) that is 
introduced in Section 4.1. To show that the limit lim^^ exists, we insert the definition 
of Aq(£) into eq. (33) . After splitting off all non-singular terms in Aq we obtain: 



Sq(z) = lim 



2 Re 



= lim 

V— >x 



( A (£) - log(u - v) - Re log P(v) 

J v 

1 



2 Re 



d£ - log(v - v ) + reg„_ 



lim [ log(i> 



log(z — z) — log(-i> — v) + reg' 



Since the singularity from the integral cancels against the term log(v — v), the limit can 
be taken. 

Our second aim is to understand that Sq(z) = Sq (z) does not depend on the choice of x. 
Let us displace x by some small amount o 6 1 such that x + a is still in the Dirichlet 
interval D\. Comparison of Sq (z) and Sq +(1 (z) gives 

S*(z)-S* + a (z) 



lim 

v— >x 



lim 



lim 

v— >x 



2 Re 



-2Re 



v-\-a 



A (£)^-RelogP(» +RelogP(v + a) 



>+a 2 9~1 



E 



+ Re log 



P(v + a) 



23+1 



-Re ^ (log(v + a - Xi) - \og{v - Xi)) + Relog 



i=l 



P(v + a) 







In passing to the second line we omitted all terms in the integrand which vanish when £ 
comes close to the real axis. 

Finally, we investigate the behaviour of So(z) at the boundary. Basically, one repeats 
the analysis we have sketched above in our discussion of lim^^^. If the end-point z of 
our integration approaches one of the Dirichlet intervals, this leads to the singularity 
~ — log \z — z\. In the argument one needs that the quotient a/ P(z) / P(z) in front of the 
singular term l/(z — z) satisfies lim z ^ x a/ P(z)/P(z) = 1 for x G Dj~. This is no longer 
true when z is sent to the real axis in one of the Neumann intervals Nk- In fact, upon 
moving x from a Neumann into a Dirichlet interval, the polynomial P(x) changes sign, 
causing the quotient P(z)/P(z) to surround the origin of the complex plane once. After 
taking the square root we conclude that lirn^^ ^JP{z)/P{z) = — 1 for x £ and hence 
Sq(x) ~ log \z — z\ near the Neumann intervals. 
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